Amplitude modulated drift wave packets in a nonuniform magnetoplasma 
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(N . 
H Abstract 

g We consider the amplitude modulation of low-fluency, lo„ g wave,e„ g ,l, electrostatic drift wave packets in a nonuniform ma B - 
_^ netoplasma with the effects of equilibrium density, electron temperature and magnetic field inhomogeneities. The dynamics of 
^ the modulated drift wave packet is governed by a nonlinear Schrodinger equation. The latter is used to study the modulational 
-} instability of a Stoke's wave train to a small longitudinal perturbation. It is shown that the drift wave packet is stable (unstable) 

qq against the modulation when the drift wave number lies in < k < 1/ y2 (1/ y2 < k < 1). Thus, the modulated drift wave packet 

t— I can propagate in the form of bright and dark envelope solitons or as a drift wave rogon. 
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' A nonuniform magnetoplasma supports a great variety 
55 of low-frequency electrostatic and electromagnetic drift -type 
modes. Examples include the electrostatic drift waves [1] and 
O ^oupled drift- Alfven waves 12], which play a crucial role in 
^ cross-field plasma particle transports U, and the formation 
Q of coherent structures in space fl and laboratory (0 0] plas- 
" 5a m as that are magnetized. Both the drift- Alfven waves can be 
>~>excited by free energy sources that are stored in the equilibrium 
i-C pressure gradient and in magnetic field inhomogeneity. Non- 



ON 



thermal drift waves attain large amplitudes and start interacting 
among themselves. In the past, Tasso 18Q and Orevskii et al. 
|@] considered nonlinear interactions between one-dimensonal 
drift waves propagating in a direction orthogonal to the den- 



sity and temperature gradients and a uniform magnetic field in 
an electron-ion plasma. They reported the formation of non- 
envelope drift solitary pulses that were used in the study of drift 
Wave turbulence comprising an ensemble of drift wave solitons 

oil 

in magnetized plasmas. 
CN| ' Hasegawa and Mima [11] incorporated the vector nonlinear- 
r j~^ aty associated with the nonlinear ion polarization drift in the 
study of nonlinearly interacting pseudo-three-dimensional drift 
k> iWaves m a nonuniform plasma without the electron tempera- 
r*jjj ture and magnetic field inhomogeities. The Haisegawa-Mima 
C3 equation, which similar to the Charney equation [ 12] governing 
the dynamics of the Rossby waves in the atmosphere, admits a 
Larichev-Rezhnik type-dipolar vortex yl [l3l LL4J, [l5|] as one of 
the possible stationary states of the drift wave turbulence. The 
mode couplings between finite amplitude drift waves can also 



generate convective cells and zonal flows 11611 . The latter pro- 
vide a better plasma confinement, since they act as a barrier for 
inhibiting the transport of the plasma particles across the exter- 
nal magnetic field direction. Tynan et al. Il7ll have presented 
an elegant review of experimental drift wave turbulence studies. 

In this Letter, we discuss the properties of modulated one- 
dimensional drift wave packets in a nonuniform magneto- 
plasma with the effects of equilibrium density, electron tem- 
perature and magnetic field gradients. It is shown that the dy- 
namics of the modulated drift wave packet is governed by a 
nonlinear Schrodinger equation (NLSE), which depicts the for- 
mation of dark and bright solitons, as well as drift rogue waves 
(or drift rogons). 

Let us consider a nonuniform electron-ion plasma in the pres- 
ence of equilibrium density and electron temperature inhomo- 
geneities in a nonuniform external magnetic field zBo(x), where 
z is the unit vector along the z- axis of a Cartesian coordinate 
system, and Bo is the strength of the magnetic field. Thus, at 
equilibrium, we have 



d_ 

dx 



n (x)T e (x) + 



Bl(x) 



8^ 



0. 



(1) 
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where no(x) and T e (x) are the unperturbed electron number den- 
sity and the electron temperature, which have gradients along 
the x-axis. 

In the presence of the low-frequency (in comparison with the 
ion gyro frequency a> c j = eBo/nijC, where e is the magnitude of 
the electron charge, m, is the ion mass, and c is the speed of 
light in vacuum), long wavelength (in comparison with the ion 
thermal gyroradius p, = Vn/^ch where Vn = {k B Tilmi) 112 is 
the ion thermal speed, ks is the Boltzmann constant, and T, is 
the ion temperature) electrostatic field (E = — V0), where <p is 
the electrostatic potential of the drift waves, the perpendicular 
(to z) component of the electron and ion fluid velocities are, 
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respectively, 



U t _L = 



c ^ „ ck B TJx) „ _ 

Z X TTT~^ 7T Z X " n el 



B (x) 

C 



eB (x)n (x) 



(2) 



Bob) 



zxVf 



— - — I — + V £ -V|V I 

B (x)co ci \dt 



(3) 



where we have assumed that \d/dt\ <K v c( - <K co ce , with d/dt = 
d/dt + Y B ■ V. Here v e , is the electron-ion collision frequency, 
co ce = eBo/m e c is the electron gyrofrequency, and m e is the 
electron mass. Furthermore, without loss of generality, we have 
taken <k T e . The electron density perturbation is denoted by 
n e \ (« n )- 

Inserting Eq. (0 into the electron continuity equation and 
using the parallel component of the inertialess electron mo- 
mentum equation with \dn( e \)/dt\ <K n el n e \, we obtain for 
v ei \dn e i/dt\ <k Vj e \d 2 n el /dz 2 \ and c[d(n /B )/dx]d 2 <p/dydt «: 
{noek B T e /v e i)\d 2 (pldz 2 \, the Boltzmann law for the electron 



number density perturbation |14 1 



n e \ = n (x) exp 

no(x) 
k B T e (x) 



k B T e (x) 
1 

d>+ - 



np(x) 
2 \kBT e (x) 



„2 jl2 



(4) 



Furthermore, substituting Eq. ((3} into the ion continuity equa- 
tion, we obtain 



dn n cn Q dV 2 ± 4> 

— + V ■ [(n (x) + n el )\ E ] * -j— 

at B()(jL> C i dt 



(5) 



where we have neglected the parallel (to z) ion dynamics, 
thereby discarded the coupling between the drift and ion-sound 
waves. 

We can now combine Eqs. © and (O under the quasi- 
neutrality condition tin = n e \, which holds for a magnetized 
plasma with oj p j » a> C i, where u> p j = {AnriQe 2 /m/) 1 ! 2 is the 
ion plasma frequency, to obtain the drift wave equation in one- 
space dimension 



d 2 \ dip dip dip dip 



(6) 



where ip = e<p/k B T e , a = p s /L n , and y - p s /L TB . Here 
Ps = C s /a) C i is the ion sound gyroradius, C s = (k B T e /mi) 1 ^ 2 
is the ion sound speed, L~ l = -dln(no/Bo)/dx > 0, and 
Lj B - dln(rio/T e Bo)/dx. Furthermore, the time and space vari- 
ables are in units of the ion gyroperiod 1 /w f , and p s . 

Let us now derive the governing nonlinear equation for 
amplitude-modulated drift wave packets, following the stan- 
dard multiple-scale technique lfl8l 19 . 20ll . Then in a coordi- 



of perturbation. We are interested in the modulation of a plane 
drift wave as the carrier wave with the wave number k and the 
frequency u>. The dynamical variable <p can be expanded as 

DO OO 

v = Z Z r) cxp[i(ky - wf)L (7) 



where ip_, = ip^"!* is the reality condition and the asterisk de- 
notes the complex conjugate. 

Substituting the expansion, given by Eq. (0, and the 
stretched coordinates in Eq. ©, and equating different pow- 
ers of e, we obtain for n — I — 1 (coefficient of e) the linear 
dispersion relation for drift waves 



(n> 



/ = — OO 



ak 



(8) 



1 +k 2 ' 

From the second-order expressions for n — 2,1 — 1 we ob- 
tain an equation in which the coefficient of ip^ vanishes due 
to the dispersion relation, and the coefficient of dtp\ /d£, after 
equating to zero, gives the group velocity 

a(l-k 2 ) 

v * = (TW (9) 

Next, the zeroth harmonic mode, which appears due to the 
nonlinear self-interaction of the carrier waves in the coefficient 
of e 3 for n = 2, 1 = 0, is obtained as 



(10) 



Considering the second order harmonic mode for n = I — 2, we 
obtain from the coefficient of e 2 



(2) 



to + yk 



(l) n 2 



(ID 



2ak-u(\ + 4k 2 ) 
Finally, for n = 3,1 = 1, we obtain an equation for third - 
order first-harmonic modes in which the coefficients of ip\ and 
dipf } jdt; vanish by the dispersion relation and the group ve- 
locity, respectively. In the reduced equation we substitute the 
expressions for ^[ 2) and ip? from Eqs. (fTUl i and (fTTT i. Thus, we 
obtain the following NLSE 



SO 5 2 4> 
' — + P 

8t df 



i^- + + Q\®\ 2 ® = 0, 



(12) 



where O = <p\ is the potential perturbation, or in terms of the 
original variables 

i (l + v 4) a> + P + QIO|2 ° = °' (13) 

where O ~ &p\ . The coefficients of the drift wave group dis- 
persion and the nonlinearity are 



nate frame moving with the speed v g , the space and the time an( j 
variables can be stretched as £ = e(y - v g t), r = e 2 t, where e 
is a small parameter (0 < e <K 1) representing the weakness 



_ 1 d 2 u) _ ak(3 - k 2 ) 
= 2'dic 2 ~ ~ (1 + k 2 f ' 

[a + y(l+k 2 )]Q 



ak{\ +k 2 ) 2 (3 + k 2 )(l -2k 2 )' 



(14) 



(15) 



where Q Q = a(l + 3k 4 ) + y [(1 + k 2 ) 2 + k 2 (l - k 4 )]. 

The propagation of wave packets in a dispersive nonlinear 
plasma medium has been known to be subjected to the ampli- 
tude modulation, i.e., a slow variation of the wave packets en- 
velope due to the nonlinear self-interaction of the earner wave 
modes. The systems evolution is then governed through the 
modulational instability (MI). The latter signifies the exponen- 
tial growth of a small plane wave perturbation as it propagates 
in plasmas. The gain leads to amplification of sidebands, which 
break up the otherwise uniform wave and lead to energy local- 
ization via the formation of localized structures. Thus, the MI 
acts as a precursor for the formation of bright envelope soli tons, 
in absence of which we have the formation of dark solitons. 

Let us now consider the amplitude modulation of a plane drift 
wave solution of Eq. dT2b of the form <£ = ®oe~ ,n ° T , where 
Qo = -Q^>q with <Do denoting the potential of the drift wave 
pump. We then modulate the drift wave amplitude as a plane 
wave perturbation with frequency Q. and wave number K as 
O = (<D + <5>ie iK t- ,at + O 2 e" iX?+iQr )e- in0T , where (D u are real 
constants. Looking for the nonzero solution of the small per- 
turbations, we obtain from Eq. ([T21 the following dispersion 
relation for the modulated drift wave packets: 



(16) 



where K c - ~\/2\Q/P\\Oo\ is the critical value of the wave num- 
ber of modulation K, such that MI sets in for K < K c , and the 
wave will be modulated for PQ > 0. In the latter, the perturba- 
tions grow exponentially during propagation of waves. On the 
other hand, for K > K c the wave is said to be stable {PQ < 0) 
against the modulation. The instability growth rate is obtained 
as 



(17) 



Clearly, the maximum value of F is achieved at K = K c / V2 
and is given by r max = l2l|Oo| 2 . From Eqs. ([Pil l and (fT3T > 
we note that for long wavelength drift modes (with k < 1), 
the dispersive coefficient P is always negative, whereas Q ^ 
according as k «s 1/ y2 » 0.71. Thus, the drift wave packet is 
stable (PQ < 0) or unstable (PQ > 0) against the modulation 
according to when < k < 0.7 1 or 0.7 1 < k < 1. 

Next, we numerically investigate the MI growth rate as 
shown in Fig. 1. For a fixed a and y (e.g., a — y = 0.05) 
we note that as k approaches 1, the value of F decreases with a 
lower cut-off at a lower wave number of modulation. However, 
for long wavelength perturbations with K < 1, the MI growth 
rate can not be controlled for drift waves with wave numbers k 
close to 1 / y2. The reason is that when k approaches 1 / y2, the 
nonlinear coefficient Q becomes larger and larger. 

Exeat solutions of the NLSE (TT~2T > can be obtained by consid- 
ering <t> = V*Pexp(/0), where *P and 6 are real functions to be 
determined (see for details, e.g., Refs. ll2lln . For PQ > the 
drift wave is modulationally unstable leading to the formation 
of bright envelope modulated wave packets given by (Fig. 2) 



*F = *F sech 2 ( - UT ) , 6 = — 
\ W ) 2P 



(18) 




Figure 1: The instability growth rate given by Eq. dl7l is shown at different 
wave numbers of the carrier drift mode. The parameter values are a = y = 0.05. 
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Figure 2: The evolution of the bright envelope soliton as given by Eq. {TD at 
T = for k = 0.71, or = y = 0.05, *P = 0.04 and U = 0.5. 



which represents a localized pulse traveling at a speed U and 
oscillating at a frequency Qo at rest. The pulse width W is re- 
lated to the constant amplitude *Fo as W — V2f / Q^o- 

On the other hand, for PQ < 0, the modulationally stable 
drift wave packet will propagate in the form of a dark envelope 
soliton characterized by a depression of the drift wave potential 
around £ = 0. This is given by (Fig. 3) 



*F = *F, tanh 



g - Ut 



1 

2P 



■ u 2 

U^-\ — -2PQ^ 



(19) 



representing a localized region of hole (void) traveling at a 
speed U. The pulse width W2 depends on the constant ampli- 
tude *Pi as W 2 = y/2\PjQ\Ni- 

Furthermore, the NLSE (flZt has a rational solution that is lo- 
cated on a non-zero background and localized both in the t and 
£ di rections. For PQ > we have the drift rogue waves/rogons 



O = cD 



A(\+HQt) 



1+4Q 2 t 2 +2QP/P 



- 1 



exp(/£r), (20) 



Figure 3: The evolution of the dark envelope soliton as given by Eq. {TgJ at 
T = for k = 0.5, a = y = 0.05, ¥ = 0.04 and U = 0.5. 



which reveals that a significant amount of energy is concen- 
trated in a relatively small area in space. The typical form of 
the rogue wave is shown in Fig. 4. Hence, a random pertur- 
bation of the drift wave amplitude will grow on account of the 
modulational instability. 

To summarize, we have considered the amplitude modulation 
of a finite amplitude one-dimensonal electrostatic drift wave 
packet in a nonunform magnetoplasma in presence of equi- 
librium density, electron temperature and magnetic field gradi- 
ents. It is shown that the dynamics of the modulated drift wave 
packet is governed by a nonlinear Schrodinger equation. Since 
the group dispersion of the drift waves is always negative for 
Ic < 1, the formation of the bright envelope drift wave soliton 
or drift rogons is possible only when Q < 0, which happens for 
0.71 < Ic < 1. In the opposite case, viz. < k < 0.71 when 
Q > 0, an amplitude modulated drift wave packet is stable and 
it propagates in the form of a dark envelope soliton |24J. In 
conclusion, the present results should be helpful in identifying 
modulated drift wave packets that may spontaneously emerge 
in magnetized space and laboratory plasmas that contain equi- 
librium density, electron temperature and magnetic field inho- 
mogeneities. 
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